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ABSTRACT 

We present X-ray /7-ray spectra of the binary GX 339-4 observed in the hard state 
^ ' simultaneously by Ginga and CGRO OSSE during an outburst in 1991 September. 

The Ginga X-ray spectra are well represented by a power law with a photon spectral 
index of T ~ 1.75 and a Compton reflection component with a fluorescent Fe Kcv line 
corresponding to a solid angle of an optically-thick, ionized, medium of ~ 0.4 x 2ir. The 
OSSE data 50 keV) require a sharp high-energy cutoff in the power-law spectrum. 
The broad-band spectra are very well modelled by repeated Compton scattering in a 
thermal plasma with an optical depth of r ~ 1 and kT ~ 50 keV. We also study the 
distance to the system and find it to be 3 kpc, ruling out earlier determinations of 
' ~ 1 kpc. Using this limit, the observed reddening and the orbital period, we find the 

allowed range of the mass of the primary is consistent with it being a black hole. 
q ■ We find the data are incosistent with models of either homogenous or patchy 

5h ' coronae above the surface of an accretion disc. Rather, they are consistent with the 

c/3 , presence of a hot inner hot disc with the viscosity parameter of a ~ 1 accreting at 

a rate close to the maximum set by advection. The hot disc is surrounded by a cold 
outer disc, which gives rise to the reflection component and a soft X-ray excess, also 
present in the data. The seed photons for Comptonization are unlikely to be due 
' to thermal synchrotron radiation. Rather, they are supplied by the outer cold disc 

J-j \ and/or cold clouds within the hot disc, pair production is negligible if electrons are 

thermal. The hot disc model, which scaled parameters are independent of the black- 
hole mass, is supported by the similarity of the spectrum of GX 339-4 to those of other 
black-hole binaries and Seyfert Is. On the other hand, their spectra in the soft 7-ray 
regime are significantly harder than those of weakly-magnetized neutron stars. Based 
on this difference, we propose that the presence of broad-band spectra corresponding 
to thermal Comptonization with kT il 50 keV represents a black-hole signature. 

Key words: accretion, accretion discs - binaries: general - gamma-rays: observations 
- gamma-rays: theory - stars: individual (GX 339-4) - X-rays: stars. 



1 INTRODUCTION 

GX 339-4, a bright and well-studied binary X-ray source, is 
commonly classified as a black hole candidate based on the 
similarity of its X-ray spectral states and short-time variabil- 
ity to those of Cyg X-l (e.g. Tanaka & Lewin 1995). How- 
ever, determinations of the mass of its compact star, Mx, 
have been inconclusive (e.g. Cowley, Crampton & Hutchings 
1987, hereafter C87; Callanan et al. 1992, hereafter C92), 



and thus its nature has been uncertain. Therefore, further 
studies of the properties of GX 339-4 as well as their com- 
parison to those of objects with more direct evidence for 
harbouring a black hole is of crucial importance. 

In this work, we present two, very similar, broad-band 
X-ray/7-ray (hereafter X7) spectra of GX 339-4 obtained 
during a strong outburst of the source in September 1991 
(Harmon et al. 1994) simultaneously by Ginga (Makino et 
al. 1987) and the Oriented Scintillation Spectroscopy Ex- 
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periment (OSSE) detector (Johnson et al. 1993) on board 
the Compton Gamma Ray Observatory (CGRO). The source 
was in the hard (also called 'low') spectral state. The Ginga 
and OSSE observations were reported separately by Ueda, 
Ebisawa & Done (1994, hereafter U94) and Grabelsky et al. 
(1995, hereafter G95), respectively. However, the data from 
the two instruments have not been fitted together, and, e.g. 
G95 found models with Compton reflection of X7 photons 
from an accretion disc unlikely whereas U94 found strong 
evidence in the Ginga data for the presence of this process. 

Here, we re-analyze the simultaneous Ginga and OSSE 
data based on the present accurate calibration of those in- 
struments. This leads to a reconciliation of the apparent dis- 
crepancies between the data sets from the two instruments, 
and allows us to fit the joint data with physical models. We 
also study the distance, reddening, Galactic column den- 
sity and the masses of the binary members. Those results 
are then used in studying radiative processes, geometry and 
physical models of the source. Finally, we find the X7 spec- 
trum of GX 339-4 similar to those of black-hole binaries 
and Seyfert AGNs, and, in particular, virtually identical to 
that of NGC 4151, the Seyfert brightest in hard X-rays. This 
favours physical models with scaled parameters independent 
of the central mass, such as a hot accretion disc with unsatu- 
rated thermal Comptonization (Shapiro, Lightman & Eard- 
ley 1976, hereafter S76). On the other hand, the spectrum of 
GX 339-4 is significantly different from those observed from 
neutron star binaries, which supports the black-hole nature 
of the compact object in GX 339-4. 



2 THE PARAMETERS OF THE BINARY 

In order to analyze the X-ray data meaningfully, we need 
to estimate basic parameters of the binary system. Of im- 
portance here are the Galactic column density, Nh, the in- 
terstellar reddening, Eb-v, the distance, d (for which pub- 
lished estimates range from 1.3 to 4 kpc), the masses of the 
primary and secondary, Mx and M c , respectively, and the 
inclination (with respect to the normal to the orbital plane), 
i. 



2.1 Reddening and column density 

Grindlay (1979) found strong interstellar Nal D absorption 
lines and diffuse interstellar bands at A ~ 5775-5795, 6010, 
6176, 6284, and 6376 A, while C87 found a strong interstel- 
lar Call K absorption line and diffuse A4430 A absorption 
band. The equivalent widths of these features are consistent 
with Eb-v — 1—1.3. From the uncertainties of the published 
estimates, we derive the weighted mean of 

Eb-v = 1-2 ± 0.1 . (1) 

The most extended all-sky study of the distribution of 
neutral H based on high-resolution IUE observations of Lya 
absorption towards 554 OB stars shows their Nh well corre- 
lated with the column density of dust, measured by Eb-v, 
with (N H / 'Eb-v) = 4.93 x 10 21 cm" 2 mag" 1 (Diplas & Sav- 
age 1994). Eb-v = 1.2 ± 0.1 derived above thus indicates 

TVh = (6.0 i 0.6) x 10 21 cm -2 . (2) 



This Nu is in excellent agreement with that derived 
from X-ray data. We obtain N H = (6.2 ± 0.7) x 10 21 cm" 2 
from the depth of the O edge of to = 2.6 ± 0.3 measured 
by Vrtilek et al. (1991), and assuming the O abundance of 
Anders & Ebihara (1982). On the other hand, Vrtilek et 
al. (1991) and Ilovaisky et al. (1986) have obtained Nh = 
(6.6 ± 0.3) x 10 21 cm -2 and (5.0 ± 0.7) x 10 21 cm -2 from 
continuum fits in the soft and hard state, respectively. Those 
values are less reliable because those authors assumed the 
continuum models of optically-thin thermal bremsstrahlung 
and a single power law for the corresponding two states. 
The first model cannot hold for a luminous source based 
on the standard efficiency argument, and the second model 
disagree with Ginga data showing the presence of a strong 
soft X-ray excess in addition to the harder power law over 
two orders of magnitude of the flux in the hard state (U94). 
Therefore, we assume hereafter JVh = 6 x 10 21 cm" 2 , which 
value agrees with both optical and X-ray data. 

On the other hand, Ilovaisky et al. (1986) obtained 
Eb-v = 0.7 ± 0.1 by converting their fitted Nh (which ap- 
pears underestimated, see above) with an old calibration 
of Nh/Eb-v, which underestimates Eb-v(Nh) according 
to the presently most extensive study of Diplas & Savage 
(1994). This value of Eb-v is also in disagreement with the 
results from interstellar features, and it appears incorrect. 

2.2 Distance 

Obviously, reddening increases with distance. However, the 
resulting correlation depends sensitively on direction be- 
cause the distribution of the interstellar matter is very com- 
plex and patchy (e.g. Diplas & Savage 1994; Neckel & Klare 
1980), especially towards the Galactic Center, in which di- 
rection is the line of sight of GX 339-4. 

To illustrate this problem, we show here examples of a 
wide dispersion of Es-v(d) with l,b within ±5° from the 
direction to GX 339-4. That field includes many OB stars 
in stellar systems, for which both Eb-v and d have been 
well studied. First, the distribution of OB stars and dust 
m a 21-deg 2 field centered at (/,&) = (335°, 0°) was stud- 
ied in detail by FitzGerald (1987), who found the dust is 
distributed in two distinct clouds, one in the local arm at 
190 ± 30 pc and the other in an interarm cloud at 690 ± 70 
pc with (Eb-v) = 0.21 and 0.76, respectively. Three OB 
associations along this line of sight found at (d) = 1.34, 2.41 
and 3.69 kpc have Eb-v = 0.82, 0.86 and 0.92, respectively. 
On the other hand, the open clusters NGC 6200, NGC 6204 
and Hogg 22 at (I, b) ~ (338.5, -1.2) show E B -v < 0.66 for 
d < 2.78, while NGC 6193 at (I, b) = (336.7, -1.6) and NGC 
6250 at (340.8,-1.8) show E B -v = 0.44-0.49 at d = 1.4 kpc, 
and Eb-v = 0.38 at d — 0.95 kpc, respectively (Moffat & 
Vogt 1973, 1975; FitzGerald et al. 1977; Vazquez & Fein- 
stein 1992). Finally, the clusters NGC 6208 at (333.5,-5.7), 
IC 4651 at (341,-8) and NGC 6352 at (342,-7) at d = 1, 
1.15 and 5 kpc show Eb-v = 0.18, 0.12 and 0.4, respectively 
(Lindoff 1972; Alcaino 1971). These results clearly show that 
EB-v(d) around GX 339-4 depends sensitively on direction. 

Therefore, we have determined Eb-v and d for ~ 450 
stars with V 10 within ±5° around GX 339-4 using 
their UBV magnitudes and spectral types from the cata- 
logue of Mermilliod & Mermilliod (1994). Our study con- 
firms the strong dependence of Es-v(d) on direction and 
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Figure 1. The dependence of Eb—v on distance in the field 
centered at GX 339—4. Dots represent stars in a small field with 
a radius of ~ 1° around the source. Two stars closest to the line 
of sight of GX 339-4 are shown by the two dots between 2 and 

3 kpc. Crosses and open circles represent nearby sources outside 
that radius on the side of the Galactic Center and on the opposite 
side, respectively. (Note an increase of extinction per unit distance 
toward the Galactic Center.) On the other hand, the extinction 
of GX 339-4 is E B -v = 1.2 ± 0.1. 

allows an accurate determination of EB-v(d,l,b) in that 
field. We have found that Eb-v depends much more sen- 
sitively on b than on I, and thus we present in Fig. 1 our 
results only for a smaller field centered at GX 339-4 with 
Q,b) = (338.9 ± 5°, -4.3 ± 1°). Stars with I <> 338° and 
<; 340° are indicated by open circles and crosses, respec- 
tively. Dots represent stars in a central field with a radius 
of ~ 1°. From Fig. 1, it is clear that the extinction in- 
creases up to (Eb-v) ~ 0.6 at d « 2 kpc. There is a 
possible further increase beyond ~ 3 kpc, but to confirm 
it we would need a deeper survey than that of Mermil- 
liod & Mermilliod (1994). Two stars nearest to GX 339- 

4 with known UBV magnitudes and spectral types are 
the early B-type giants HD 153222 and CPD-48°11361 at 
(/, b) w (338°, -4°), for which we have found E B -v = 0.53, 
d <; 2.2 kpc, and Eb-v = 0.54, d ~ 2.6 kpc, respectively. 
Thus, our strong conclusion is that GX 339-4 with the mea- 
sured Eb-v = 1.2 ± 0.1 must be located at d > 2.6 kpc. 

This result rules out the most recent estimate of d = 1.3 
kpc towards GX 339-4, obtained from a ROSAT measure- 
ment of its X-ray halo by Predehl et al. (1991). They suggest 
that the halo is produced by a dust cloud near the source, 
and identify that cloud with the second of two dense clouds 
found at d ~ 100-250 pc and ~ 800-1400 pc, respectively, 
by Neckel & Klare (1980). However, we have not been able to 
reproduce their identification. There are indeed two clouds 
shown in the direction of / = 339° in Fig. 9 of Neckel & 
Klare (1980), but that figure gives the cloud distribution pro- 
jected onto the Galactic plane. In fact, Table IVa of Neckel 
& Klare (1980) shows that the field of GX 339-4 (numbered 
196 in that paper) does not contribute to extinction in those 
clouds, as well as Fig. 8c of that paper shows that Av at 1 
kpc in that field is < 1.9 (whereas Av ~ 3.4-4 corresponds 
to Eb-v = 1.2 ± 0.1 of GX 339-4). These findings are con- 
firmed by a more extensive study of Galactic extinction by 
FitzGerald (1987). She shows the same clouds are located 



in the field of Ara OB 1 association at (l,b) ~ (335°, 0°), 
confirming that they do not concide with the line of sight to 
GX 339-4. Moreover, the total E B -v ~ 0.8 at ~ 1.34 kpc 
(FitzGerald 1987) produced by the two clouds is lower than 
the Eb-v — 1.2 measured for GX 339-4 (and adopted by 
Predehl et al. 1991 in their analysis). Similarly, d = 1.33 kpc 
derived by Mauche & Gorenstein (1986) for GX 339-4 from 
Einstein observation of the X-ray halo corresponds just to 
the Galaxy's dust layer, and thus it represents a lower limit 
to the distance. That limit can be actually higher id ^ 2 
kpc) since recent studies (Diplas & Savage 1994) show a 
value of the scaleheight for the dust layer of 150 pc that is 
larger than 100 pc adopted by Mauche & Gorenstein. Those 
two scaleheights predict Eb-v ~ 0.5 and ~ 0.8, respectively, 
both lower than the E B -v of GX 339-4. 

On the other hand, the distance to GX 339-4 can be es- 
timated from its systemic velocity of Vo = —62 ± 10 km s 
(C92). Using the Vo-d conversion chart from Burton (1992), 
the kinematic distance is d — 4 ± 1 kpc. (Note that a pos- 
sible peculiar velocity of the system as well as a systematic 
uncertainty of the above Vo makes this estimate less secure 
than the extinction limit.) 

We note here that both the studies of distribution of 
matter in the Galaxy as well as the measurements of radial 
velocities of Hn regions in the range I w 300°-340° indi- 
cate the presence of a spiral arm at a distance of ~ 4 kpc 
(Scheffler & Elsasser 1987 and references therein). Our de- 
rived kinematic distance and the relatively large Eb-v of 
GX 339-4 are consistent with its location in that arm. 

Finally, we can put an upper limit on d from the Ed- 
dington limit, which appears to be satisfied in all known 
black-hole binaries (Tanaka & Lewin 1995). The most lu- 
minous state of GX 339-4 reported as yet appears to be a 
soft state in 1988 (Miyamoto et al. 1991), for which we esti- 
mate the bolometric flux (using the disc blackbody model, 
see Section 5.2) to be ~ 4 x 10 -8 erg cm -2 s _1 . This cor- 
responds to d ~ 10(L/L e ) 1/2 (M x /3Mq) 1/2 kpc, where L E 
is the Eddington luminosity (see Section 2.3 for discussion 
of Mx). On the other hand, typical maximum luminosities 
of black-hole binaries are at least factor of a few below Le 
(Tanaka & Lewin 1995), and, e.g. d ~ 4.5 kpc if L — 0.2Le 
and M x = 3M Q . 

Concluding this section, our strongest limit is d <; 3 
kpc based on the extinction study. This limit is consistent 
with the kinematic distance of d = 4 ± 1 kpc and with the 
Eddington limit. Thus, we adopt d — 4 kpc hereafter. 

2.3 The masses and geometry 

With the above results on the Eb-v and d, we can con- 
strain the mass of the companion star. In general, we expect 
M c ^ IMq in low-mass X-ray binaries (hereafter LMXBs, 
e.g. van Paradijs & McClintock 1995), to which class GX 
339-4 most likely belongs. At d ~ 4 kpc and A v ~ 3.4-4, 
a lM© main sequence star has an observed V ~ 20.5-21 
(B - V ~ 1.8, V - Rc ^ 1, V - Ic =± 2), which corresponds 
to fluxes below the faintest ones observed as yet from GX 
339-4, as illustrated in Fig. 2. Thus, the presence of a IMq 
secondary star is entirely consistent with the present data. 
We note that our conclusion differs from that of C92, who 
claimed M c ^ 0.4Mq. The reason for the discrepancy is 
that those authors assumed d = 1.3 kpc (from Predehl et al. 
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Figure 2. Fluxes at the / c , R c , V, and B bands (from left 
to right) of a IMq main sequence star at 4 kpc reddened by 
Eb — v = 1-2 (open squares) and by Eb—v = 1 (a less likely 
value; open circles) compared to the faintest observations (not 
simultaneous) of GX 339-4 (filled circles). The B, V are from Ilo- 
vaisky & Chevalier (1981) and Ilovaisky et al. (1986), and -R c , Ic 
are from C92. For both estimates of Eb—Vi the stellar spectrum 
is below the observations of GX 339—4 (which spectrum is a sum 
of emission from the companion star and from the accretion flow 
onto the compact object). Thus, the presence of a IMq secondary 
star is fully compatible with the present data. 



1991) and E B -v = 0.7 (from Ilovaisky et al. 1986), which 
both values have been found by us to be incorrect, see above. 

C92 detected a 14.8-h flux modulation in both optically- 
high and low states. They also found the same periodicity 
in the radial velocities of He n 4686 and H I Balmer emission 
lines measured in low resolution spectra from 1986 May 7-9 
(C87). We caution here that interpreting this periodicity as 
due to the orbital period remains relatively uncertain. 

First, the periodicity is absent in the radial velocity 
data from 1985 March 14-17 of C87, although it is possibly 
due to their poorer quality. Second, although C87 claimed 
that the source was in its 'normal' (high optically) state of 
~ 16.5 mag both in 1985 March and 1986 May, C92 noted a 
displacement in the mean velocity between those data from 
~ +96 km s" 1 to —62 ± 10 km s~ , as well as Ilovaisky et 
al. (1986) reported a relatively faint optical state from late 
1985 March (B ~ 18.8) through 1985 April 29, which con- 
tradict this claim of C87. Third, Corbet et al. (1987) shows 
strong variability around 1986 May. Namely, their optical 
photometry from 1986 April and June/July shows variabil- 
ity by several tenths of magnitude during each observing 
night and an overall brightening by 0.5 mag from April to 
June. On the other hand, the equivalent widths of both He n 
4686 and H/3 lines for both data sets are practically the same 
(C87), and similar to the mean equivalent widths observed 



in LMXBs (Smale 1996). However, these lines can be formed 
in either the heated hemisphere of the secondary, the accre- 
tion disc, and/or the accretion stream, depending both upon 
the object and/or its state of activity, which precludes an 
unambiguous interpretation of the spectroscopic periodicity. 

Furthermore, we note that C92 do not state whether 
the 14.8-h periodicity found by them in the (optically) high- 
state photometry of Corbet et al. (1987) was present in both 
1986 April and June/July or in only one of those data sets. 
C92 also claim that if the maximum in the light curve (pre- 
sented in their Fig. 6) corresponds to transit of the compact 
object in front of its companion, the radial velocity changes 
(shown in their Fig. 7) are consistent with an origin near 
the compact object. However, this is not the case because 
at both spectroscopic conjunctions (i.e., the transits of the 
compact star both in front and behind the companion) the 
observed radial velocity should be just the systemic velocity. 

Keeping these doubts in mind let us tentatively assume 
that the 14.8-h modulation is still due to the orbital mo- 
tion. The semi- amplitude, K x = 78 ± 13 km s" 1 (C92; the 
uncertainty is 1-a) implies the secondary mass function of 



M c 



(1 + M x /M c ) 



= 0.030] 



0.012 M e i 



(3) 



which corresponds to the mass of the primary of 

Mx = M c [(M c //) 1/2 sin 3/2 i - l] . (4) 

From the lack of eclipses, the inclination can be con- 
strained to t £ 60° (C87). If we adopt M c = 1M Q (see 
above), we obtain Mx ^ 5Mq, where the equality corre- 
sponds to the minimum / and maximum i. On the other 
hand, the lack of X-ray dips and fits of models with Comp- 
ton reflection (see Section 4) both favour < 60°, and, e.g., 
i — 45° yields Mx 3.4Mq. Current theoretical (see, 
e.g. Haensel 1995 for a review) and observational (e.g. van 
Paradijs & McClintock 1995) estimates yield the maximum 
mass of a neutron star of ^ 2Mq. On the other hand, if we 
assume that the compact object is a black hole, as strongly 
suggested by the similarity of its spectral and timing be- 
haviour to Cyg X-l, then the presumed Mx <; 2Mq corre- 
sponds to M c <; O.5M0 (at the lower limit of /). Thus, the 
mass function of C92 together with the constraints on M c are 
fully consistent with the presence of a black hole in GX 339- 
4 (although the presence of a neutron star is also allowed). 
Hereafter, we will adopt fiducial values of Afx = 3Mq and 
i — 45° , which are around the middle of the allowed param- 
eter space. 

Then, P orb = 14.8 h and M x + M c ~ 4M imply a 
1.3-Rq Roche-lobe (tidal) radius for the secondary. A post- 
main sequence IMq star with R ~ 1.3-Rq (and thus fill- 
ing its Roche lobe) would have a luminosity of ~ ILq and 
an effective temperature (unheated) of ~ 5100 K (Web- 
bink, Rappaport & Savonije 1983). Such star at d — 4 
kpc and E B -v = 1.2 would have V ~ 21.6, B - V ~ 2, 
V — i?c ~ 1.1, and V — I c ~ 2.3, which correspond to 
fluxes below those observed from GX 339-4 (similarly to 
the case of a IMq main-sequence star, see Fig. 2). The ex- 
pected Roche-lobe overflow rate, M, in such a system is 
several times 10" 10 M Q yr" 1 (Webbink et al. 1983), which is 
sufficient to power the average X7 emission from GX 339-4 
(Rubin et al. 1998) of ~ 2 x 10 36 erg s _1 assuming an ac- 
cretion efficiency of r\ = 0.06. On the other hand, M can be 
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significantly higher due to X7 irradiation of the companion 
star (Podsiadlowski 1991; Tavani & London 1993), which 
would imply a lower r\. 

In this geometry, the secondary star subtends a solid an- 
gle of ~ 0.018 x47r. Thus, the secondary provides a negligible 
contribution to Compton reflection in GX 339-4 (Section 4). 



3 Ginga AND OSSE DATA 

GX 339-4 was observed by Ginga on 1991 September 11 
and 12 (U94), and by OSSE 1991 September 5-12 (G95). 
We have extracted the OSSE data for September 11 and 12 
for periods approximately coinciding with the Ginga obser- 
vations. The log of the Ginga (from U94) and OSSE obser- 
vations used here is given in Table 1. Both the Ginga and 
OSSE fluxes during September 11-12 were approximately 
constant (U94; G95). 

We adopt the current best estimate of the effective area 
of the Large Area Counter (LAC) of Ginga (Turner et al. 
1989) of 4661 cm 2 (D. Smith, private communication). The 
usable energy range of Ginga is 1.2-29 keV for these obser- 
vations. A 1 per cent systematic error is added in quadrature 
to the statistical error in each Ginga channel (as in U94). 

The OSSE data are from 50 keV to 1000 keV. They 
include energy-dependent systematic errors estimated from 
the uncertainties in the low-energy calibration and response 
of the detectors using both in-orbit and prelaunch calibra- 
tion data, and correspond to an uncertainty in the effective 
area in the OSSE response. They are most important at the 
lowest energies (~ 3 per cent at 50 keV, decreasing to ~ 0.3 
per cent at J> 150 keV). 



4 ELEMENTARY SPECTRAL MODELS 

In this section, we analyze the Ginga and OSSE data us- 
ing idealized models of blackbody emission, and power-law 
or Comptonization emission from a hot, isotropic, thermal 
plasma cloud including their Compton reflection from an un- 
derlying slab. Here, we assume the reflected component is 
not Comptonized by the hot plasma. This allows us to char- 
acterize spectral components present in the X7 spectrum in 
a phenomenological, but relatively model-independent, way. 
Similar models are widely used in literature to model spec- 
tra of other black- hole sources, e.g. Cyg X-l (Ebisawa et al. 
1996) or Seyfert Is (Nandra & Pounds 1994), and thus the 
results of this Section allow a direct comparison of the ob- 
tained parameters with those for other objects. On the other 
hand, we introduce more realistic geometries in context of 
specific models and treat Comptonization of a fraction of 
the refl ected spectrum as well as energy balance in Sections 



5.2 below. 



For spectral fits, we use xspec (Arnaud 1996) vlO. The 
confidence ranges of each model parameter are given for a 
90 per cent confidence interval, i.e., A\ 2 = 2.7 (e.g. Press 
et al. 1992). On the other hand, the plotted vertical error 
bars are 1-cr, the upper limits, 2-a, and the plotted spectral 
data are rebinned for clarity of the display. Model spectra 
are attenuated by JVh = 6 x 10 21 cm -2 (Section 2.1). Since 
the Ginga exposure is much longer for the data set 1, we 



discuss below results obtained with that set. However, we 
give fit results for both data sets in Table 2. 

As found by U94, the Ginga spectra of GX 339-4 in the 
hard state consist of four main components: an underlying 
power law, a soft excess below ~ 4 keV, a continuum due 
to reflection from the surface of an ionized accretion disc 
(e.g. Lightman & White 1988; George & Fabian 1991), and 
a fluorescent Fe Ka line. 

For modeling Compton reflection, we use inclination- 
dependent Green's functions of Magdziarz & Zdziarski 
(1995) (instead of the angle-averaged reflection spectrum 
used in U94), which assumes an isotropic point source (or, 
equivalently, an optically-thin corona) above a slab. The 
Green's functions are convolved with an incident continuum 
(a power law or a thermal Comptonization spectrum). The 
reflector inclination is kept at i — 45° (unless stated other- 
wise). We treat the solid angle, fl, subtended by the reflector 
as a free parameter. Values of Q < 2n may correspond ei- 
ther to a truncation of the reflecting disc, or to substantial 
Comptonization of the reflected spectrum (or both). The 
ionization parameter of the reflector, £ = Li on /nr 2 , is as- 
sumed to be uniform in the reflecting material. Here Lion is 
defined as the 5 eV-20 keV luminosity in a power law spec- 
trum and n is the density of the reflector located at distance 
r from the illuminating source (Done et al. 1992). The re- 
flector temperature is kept at 10 6 K, which is the highest 
temperature consistent with the model of ionization equilib- 
rium used (Done et al. 1992) . This temperature is consistent 
with our estimate of the origin of reflection from a larger area 
than that giving rise to the observed soft excess, see Section 
5.4. The abundances of both the reflector and the interstellar 
medium are from Anders & Ebihara (1982) except that the 
relative Fe abundance in the reflector, Af b , is a free parame- 
ter. The ion edge energies and opacities are from Reilman & 
Manson (1979, used in U94) except that now the Fe K-edge 
energies are from Kaastra & Mewe (1993). The continuum 
reflection is accompanied by an Fe Ka line, which we model 
as a Gaussian centered at an energy, Ef b , with a width, ap c , 
and a flux, 7f . 

We first fit the data set 1 from Ginga in the 4.1-29 keV 
range only, where a contribution from the soft excess was 
found to be negligible by U94. We first use a model consist- 
ing of a power law and a line. This model provides a very 
poor description of the data, with x 2 = 209/30 d.o.f. The 
pattern of residuals is characteristic of the presence of Comp- 
ton reflection, as shown in Fig. 2a of U94. Indeed, adding a 
reflection component improves the fit dramatically, reducing 
X 2 to 11.7/27 d.o.f., see fit la in Table 2 (which also give cor- 
responding results for the data set 2) and Fig. 3. Compton 
reflection is present at very high statistical significance, with 
a probability of < 10 _lD that adding this component were 
not required by the data (as obtained using the F-test). The 
statistical significances of the reflector being ionized, over- 
abudant in Fe, and an Fe Ka being present in the spectrum 
correspond to the probability that the resulting fit improve- 
ment were by chance of < 10 -4 , 3 x 10~ 7 and 3 x 10 -4 , re- 
spectively. The width of the line is only weakly constrained, 
o"Fo 0.5 keV at la, and thus we keep it fixed at 0.1 keV 
hereafter. 

We comment here on the low reduced x 2 obtained for 
the data set 1. This is due to the 1 per cent systematic er- 
ror added to the statistical error. This value appears to be 
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Table 1. The log of observations in UT. The counts of Ginga and OSSE are for 1.2-29 keV and 50—150 keV energy ranges, respectively. 
The OSSE exposure and count rate are normalized to one of its 4 detectors operating during the observations. 











Ginga 








OSSE 




Data set 


Date 


Start 


End 


Livetime [s] 


Counts [s 1 ] 


Start 


End 


Exposure [s] 


Counts [s 1 ] 


1 


1991 Sept. 11 




Q5 h 26 m 


1952 


1716 ±4 


03' l 25 m 


05 h 47 m 


4795 


13.14 ±0.26 


2 


1991 Sept. 12 




01 h ll m 


384 


1695 ± 5 


01 h 10 m 


01 h 59 m 


2488 


12.71 ± 0.35 



Table 2. The parameters of the fits in Section 4 to the Ginga and OSSE data. For each data set, the consecutive fits are to the 4.1—29 
keV, 4.1—1000 keV, and 1.2—1000 keV ranges, respectively. L and L^b are the unabsorbed bolometric luminosities of the total spectrum 



and of the blackbody component only, respectively, in units of 10 37 erg s 1 . A is the 1-keV normalization of the primary continuum (in 



cm 2 s" 


- 1 keV" 


£ is in 


erg cm s 


1 r 

> J Fc 


is in 10 3 cm 


2 s 1 , kT, kT hh are in keV, and Wp e is 


in eV. 








No. 


A 


r 


kT 


L 


f2/27r 


c 

1991 


September 


kT hh 

11 


Lbb 


E Fc 






X 2 /dof 


la 
lb 
Ic 


0.75 
0.75 
0.60 


i 7C+0.02 
1 -'°-0.03 

1 76+ 02 
1 -'°-o.oi 

-, 77 +0.01 
' -0.01 




3.12 


n , 7 +0.06 
U,a '-0.05 
43+0 03 

44+O O6 
u '^-0.06 


120+ 160 
110+ 120 

OQ+90 

1991 


2 5+ 12 
z -°-0.8 

2 9 +0 1 

3 +0 
September 


25+ ' 02 
u - zo -0.03 

12 


0.29 


6 51 +°' 31 
D - Oi -0.32 

6 51+ 016 
6 54+ - 33 


1 5+°- 7 
L -°-0.8 

1 n+0.8 

1 -°-0.7 
i O+0.7 

1 -°~o.a 


49tr 6 


11.7/27 
45.2/79 
47.8/82 


2a 
2b 
2c 


0.73 
0.75 
0.55 


1 74 +0.03 
1 76+ - 02 

1 -'°-0.02 

1 76+ - 02 

1 -'°-0.03 


53 +8 7 
52tl 


2.93 


25+ 08 
u - zo -0.08 

29+ - 09 
u,zy -0.08 

29+ - 08 
u - zy -0.08 


2 o +5i0 
zou -180 

1,u -130 
170 +300 
±,u -120 


1 6+ 22 
2.0+J 
2.0±l° 


n r, 7 +0.03 
u - z -0.03 


0.21 


6 56+ ' 34 
6 "S7+0- 37 

6 5R+0- 35 
D - OO -0.38 


I 4+0-9 

i q+0.9 
L -°-0.9 
I 4+0-8 


44+Sl 

46+ 30 
TO -31 


22.4/27 
61.5/79 
62.5/82 



an overestimate of the residual inaccuracy of the calibra- 
tion of the Ginga LAC, for which 0.5 per cent appears to 
be a better estimate. However, we have retained here the 
statistical error adopted by U94 in order to enable direct 
comparison with their results. This leads to uncertainties 
on the parameters of our fits being more conservative than 
necessary, but it affects neither the best-fit values nor the 
conlusions regarding the spectral components present in the 
Ginga spectrum. E.g. Compton reflection is still found at a 
very high significance, see above. The addition of the large 
systematic error also leads, in some cases, to a relatively 
low reduced \ 2 f° r models showing systematic (rather than 
random) departures from the Ginga data, which cases we 
discuss below individually. On the other hand, the reduced 
X 2 is significantly larger for the data set 2, see Table 2. This 
is due to its much shorter exposure, resulting in typical sta- 
tistical errors > 1 per cent, which then reduces the relative 
importance of the systematic errors. 

We have then tested whether there is any indication 
in the data for the presence of kinematic and relativistic 
smearing of the Fe Ka and the reflection continuum. Such 
smearing would be present if those spectral components 
were formed close to the central black hole, which effect 
has been found in Ginga observations of the X-ray novae 
GS 2023+338 and Nova Muscae (Zycki, Done & Smith 1997, 
1998). We have used a model in which the line and the static 
reflection continuum were convolved with the disc line pro- 
file of Fabian et al. (1989), similar to the model of Zycki et 
al. (1997, 1978). However, we have found no fit improvement 
(A^ 2 = —0.6) allowing for the smearing, although the data 
also do not allow us to rule it out. 

The relative normalization of the reflected components, 
Q/2tt, corresponds to ~ 0.6 of those obtained by U94. This 
is explained mostly by their use of the angle- averaged spec- 
trum, which underestimates the actual reflection spectrum 
for angles < 65°, which effect increases at energies <; 15 
keV (Magdziarz & Zdziarski 1995). Accordingly, the fit to 



the data set 1 with the angle-averaged reflection in U94 
shows strong positive residuals above ~ 20 keV (Fig. 2b 
in U94), which systematic residuals disappear completely 
in the present fit, see Fig. 3. Also, we obtain £ about 2-3 
times less than those of U94. However, the ionization state 
depends on both £ and the reflector temperature, and the 
difference is due to the assumption of U94 that the latter 
is 10 K (which seems unlikely in inner regions of luminous 
compact objects in binaries). 

We then investigate the issue of the Fe abundance. We 
have fitted the earlier four hard-state observations of GX 
339-4 by Ginga in 1989-90 (U94), and we found all of them 
consistent with 2 ^ A Fc ^ 3, which is also the case for 
the present data. Therefore, we constrain A Fc to this range 
hereafter. 

In agreement with U94, we find a strong excess below 
4 keV, with the fluxes in the 1.2-1.7 keV and 1.7-2.3 keV 
bands about 40 and 20 per cent, respectively, higher than the 
extrapolation of the model fitted to the 4.1-29 keV range, 
as shown in Fig. 3. Although the exact form of the excess 
depends on the adopted Ah, we find it to be significant, with 
the respective 1.2-1.7 keV and 1.7-2.3 keV relative excess 
of 25 and 12 per cent even for Ah — 4x 10 21 cm -2 , which is 
much below the optical and X-ray estimates (Section 2.1). 
Thus, we conclude that the presence of the soft X-ray excess 
is not an artefact of our choice of Ah. 

At 50 keV, the extrapolated Ginga fit predict the flux 
about 10 per cent higher than that observed by OSSE. 
Above 50 keV, the OSSE data show a strong cutoff, as illus- 
trated in Fig. 3. We test whether the cutoff can be modelled 
by thermal Comptonization, similarly to the case of Cyg X- 
1 (Gierlinski et al. 1997, hereafter G97) and Seyfert AGNs 
(e.g. Zdziarski et al. 1997, hereafter Z97). We fit first the 
Ginga data above 4.1 keV jointly with the OSSE data. The 
incident continuum is fitted by thermal Comptonization in 
a spherical cloud with a central source of soft seed photons 
having a blackbody distribution, see Zdziarski, Johnson & 
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Figure 3. The ratio of the data of 1991 September 11 from Ginga 
and OSSE to the model consisting of a power law, Compton re- 
flection, and an Fe Ka line fitted to the Ginga data in the 4.1-29 
keV energy range. 



Magdziarz (1996). As discussed in that paper, the solution of 
the Kompaneets equation with relativistic corrections used 
there leads to an overestimation of the actual plasma tem- 
perature, T, when the Thomson depth of the plasma, t, is 
^ 2. Therefore, we correct here the values of T obtained from 
the Kompaneets equation by a function (R. Misra, private 
communication) obtained by comparison with correspond- 
ing Monte Carlo results (Zdziarski et al. 1996). The second 
parameter of the model is the asymptotic power-law index 
in X-rays, T, which is related to the (geometry-dependent) 
t by 



e 



-1/2 



r + ■ 



-, -1/2 



(•») 



where G = kT/m e c 2 and m c is the electron mass. In fits, we 
initially assume the seed photon temperature of fcTbb = 0.1 
keV; as long as fcTbb is much less than the minimum energy 
in the fitted data, its choice does not affect the fit results. 
This model provides a very good description of the data 
(see fits lb and 2b in Table 2). The best fit corresponds to 
r ~ 1.8. 

As an independent check of our results, we have also 
used a model of Coppi (1992), which treats thermal Comp- 
tonization using the formalism of escape probability. That 
model yields kT ~ 48 keV, r ~ 1.93, n/2-K = 0.42 as the 
best fit (x 2 = 49/79 d.o.f.), rather similar to the fit lb in 
Table 2. On the other hand, we find that the high-energy cut- 
off seen it the data is poorly modelled by an e-folded power 
law, which, apart from its mathematical simplicity, does not 
correspond to any physical model. We find that model fits 
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Figure 4. The X7 spectrum (crosses) of GX 339—4 observed si- 
multaneously by Ginga and OSSE on 1991 September 11 (the 
spectrum of September 12 is very similar and thus it is not shown). 
The data are fitted by a model consisting of blackbody radiation 
(dot-dashed curve) providing seeds for thermal Comptonization 
in a hot plasma (dashed curve). The hard radiation of the hot 
plasma is Compton-reflected from some cold medium, which com- 
ponent is shown by the dotted curve. The solid curve gives the 
sum. Note that the observed spectrum is attenuated at soft X-rays 
in the interstellar medium. The bottom panel gives the contribu- 
tion to the total x 2 from separate data bins multiplied by the 
sign of (data — model). 



the data much worse than thermal Comptonization, with 
A\ 2 = +19 resulting in a systematic pattern of residuals in 
the OSSE data. This argues for thermal Comptonization be- 
ing indeed the process giving rise to the intrinsic spectrum. 

We have considered the effect of possible inaccuracy of 
the relative normalization of the effective area of Ginga and 



OSSE. We find, however, no fit improvement, Ax 



-0.1, 



when the relative normalization of OSSE with respect to 
Ginga is allowed to be free (reaching 0.98 for both data 
sets). Therefore, we keep it at unity hereafter. 

We then check the effect of changing the disc inclina- 
tion. We find that allowing a free i leads to a negligible fit 
improvement: \ 2 = 44.3/78 d.o.f. at i = 23^23 de S- Thus, 
we keep it at i = 45° (see Section 2.3) hereafter. At the 
largest i allowed by the fit, il/2n = 0.62 and Wfc = 35 eV. 

We then consider the soft X-ray excess. In the hard state 
of Cyg X-l, the soft excess contains a blackbody component 
with fcT bb « 0.14 keV (e.g. Ebisawa et al. 1996). In addi- 
tion, the ASCA data for Cyg X-l show also a break around 
3-4 keV with the power law below the break being softer 
by AT ~ 0.4 on average (Ebisawa et al. 1996). The physical 
origin of that spectral break is unclear. We find that the 
soft X-ray excess in our data is much better fitted by an ad- 
ditional blackbody component than by a broken power law 
(Ax 2 = 5). (Our data extend down to 1.2 keV only with a 
low energy resolution, and thus do not allow us to determine 
the presence of more than one spectral component of the soft 
excess.) In the model with a blackbody component, we set 
its temperature equal to the blackbody temperature of the 
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seed photons in the Comptonization source. (We postpone 
examining the energy balance of the source to Sections 5.1- 
5.2 below.) This model gives an excellent description of the 
Ginga/ OSSF, data in the entire energy range, 1.2-1000 keV, 
see Fig. 4 and fits lc and 2c in Table 2. Still, according to 
the above discussion, the fitted value of kTbb is relatively 
uncertain; if there were a low-energy break in the power-law 
spectral component (as in Cyg X-l), kTbb would be lower. 



a) 



b) 



5 PHYSICAL IMPLICATIONS 

5.1 Geometry and energy balance 

We now consider thermal Comptonization, Compton reflec- 
tion and reprocessing in realistic geometries approximating 
the structure of an accretion flow. We take into account 
anisotropy of both seed and Compton-scattered photons, 
Comptonization of the reflected photons that return to the 
hot plasma, and energy balance. 

In general, thermal Comptonization takes place in a 
hot plasma cloud, and the resulting spectrum is Compton- 
reflected from a cold medium. A major part of the Comp- 
tonized flux incident on the cold medium is bound-free ab- 
sorbed rather than reflected, which gives rise a flux of repro- 
cessed soft photons with a spectrum close to a blackbody. 
In turn, a geometry-dependent fraction of the reflected and 
the blackbody emissions returns to the hot plasma. These 
blackbody photons then serve as seeds for thermal-Compton 
upscattering (in addition to any other soft photons present) . 
Those reflected photons that return to the hot plasma are 
upscattered, which leads to a strong suppression of that com- 
ponent for t ^ 1. 

A geometry can be ruled out either if the resulting 
model spectrum does not provide a good fit to the observed 
spectrum, or if energy balance between the hot plasma and 
the source of seed photons cannot be achieved. The condi- 
tion of energy balance can be expressed as, e.g. the power 
in soft photons incident on the hot plasma (hereafter seed 
photons) times the Comptonization amplification factor (a 
function of kT and r) being equal the power emitted by the 
hot plasma, L^ot ■ The flux in seed photons arises due to both 
internal dissipation and reprocessing of the Comptonized ra- 
diation. Then, a sufficient condition to reject a model is (i) 
the Comptonized flux (at fitted kT and r) larger than that 
observed, due to a strong flux in seed photons arising from 
reprocessing. The condition (i) is equivalent to either state- 
ment: (ii) at the fitted r and the observed hard flux, Comp- 
ton cooling by the seed photons would result in a kT less 
than that obtained from fitting, or (iii) at the fitted kT and 
the observed hard flux, energy balance can be achieved for 
values of r less than that fitted to the data (e.g. Haardt & 
Maraschi 1993; Stern et al. 1995). Below, we first obtain kT 
and t by spectral fitting, and then a posteriori check the 
energy balance by finding the equilibrium temperature cor- 
responding to the fitted r in given geometry [condition (ii) 
above] . 

In this section, we neglect spectral constraints from the 



c) 



soft X-ray excess (but return to this issue in Section 5.2 
below) because the origin of the observed soft excess is rela- 
tively uncertain, see the end of Section 4. Therefore, we con- 
sider now only the data above 4 keV, and fix the seed-photon 



d) 



e) 



Figure 5. Geometries of the X7 source considered in Sections 
5.1—5.2. Dark and light shades denote cold medium and hot 
plasma, respectively, (a) A homogeneous corona above a disc, 
(b) active regions located on the disc, (c) active regions detached 
from the disc, (d) cold clouds in the midplane of a hot disc, and 
(e) a hot central sphere surrounded by a cold disc, which may 
penetrate some distance within the sphere. Our preferred models 
are d, e. 



temperature at kTbb = 0.25 keV (as obtained in Table 2), 
which is about the highest seed-photon temperature allowed 
by the data regardless of the actual form of the soft excess. 
For a lower kTbb at a given seed-photon flux from reprocess- 
ing, the predicted Comptonized flux would increase, which 
would in turn lead to rejection of a larger class of models 
[see condition (i) above]. Thus, setting kTbb = 0.25 keV is 
the most conservative assumption for model rejection. 

To model Comptonization in anisotropic geometries, we 
use the iterative scattering method of Poutanen & Svensson 
(1996). For Compton reflection, we use the same method as 
in Section 4 (Magdziarz & Zdziarski 1995). Fig. 5 shows the 
geometries considered in Sections 5.1-{0 . 

We first consider a homogeneous corona (with the ver- 
tical optical depth of r) covering a cold slab (Haardt & 
Maraschi 1993), see Fig. 5a. In this model, all Compton- 
reflected radiation emitted by the cold slab is Comptonized 
in the hot corona, which leads to a suppression of this 
component to a level much less than that present in the 
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data. Still, the model can provide a good spectral fit if 
we allow an additional Compton-reflection component (with 
Q/2n ~ 0.3) due to an outside cold medium, e.g. an outer 
cold disc. The fitted plasma parameters are then kT = 55^* 
keV and r = 1.2±°j (x 2 = 45/79 d.o.f.). However, this 
model is strongly ruled out by the requirement of energy 
balance. Namely, the reprocessed flux from the slab is so 
strong that it would cool down the corona to ~ 27 keV. Al- 
lowing for internal dissipation in the cold slab (presumably 
a cold accretion disc) would worsen the discrepancy. 

We then consider a patchy corona geometry (Galeev, 
Rosner & Vaiana 1979; Haardt, Maraschi & Ghisellini 1994; 
Stern et al. 1995). The hot plasma is postulated here to form 
active regions above a cold accretion disc. First, we assume 
the active regions to form hemispheres located directly on 
the surface of the cold disc, see Fig. 5b. We obtain a good 
spectral fit ( X 2 = 47/79 d.o.f.) for kT = 64±g keV and 
the radial r = 2.lto'3- However, similarly to the case of a 
homogeneous corona, we find cooling by the seed photons 
arising from reprocessing in an underlying part of disc is so 
large that the plasma with the fitted r would cool to kT = 45 
keV, i.e., it cannot sustain the temperature implied by the 
data in this model. 

The cooling, however, can be reduced if the active re- 
gions are at some distance above the disc, as shown in Fig. 
5c. We find that at the height of about one radius, there 
is an energy balance between the disc without internal dis- 
sipation and the hot region. The height has to be larger if 
there is internal dissipation in the disc. However, this model 
implies much more Compton reflection (from disc regions 
surrounding the active region) than observed. Even if we as- 
sume that the underlying disc is truncated, which would re- 
duce the amount of reflection, the spectral fit of this model 
is much poorer than other models, \ 2 = 53/79 d.o.f. (for 



kT = 62±f keV, r = 2.2 



+0.2 



n/2-K = 0.71±g;$). The cause 
of the bad fit is the presence of an anisotropy break in the 
model spectrum, i.e., a defficiency of low-energy photons 
emitted upward due to an anisotropic suppression of the 
first order of Compton scattering (Haardt & Maraschi 1993; 
Poutanen & Svensson 1996). This effect is not seen in the 
data, and thus this model yields systematic departures from 
the Gmga data at low energies. On the other hand, a good fit 
can be obtained if fcTbb ^> 0.1 keV is assumed (which yields 
kT = 57±l keV, r = 2.5±gj, SI/2-k = 0.67±j};[£, X 2 = 45/79 
d.o.f.), in which case the effect of the anisotropy break on 
the model spectrum does not extend above 4 keV. The weak 
Compton reflection required by the model could then cor- 
respond to the cold disc being truncated around the active 
region, and the origin of the soft X-ray excess has to be 
due to an effect different than the disc blackbody emission. 
Thus, in this respect, the model is in principle possible. 

However, U94 found that Q. in GX 339-4 correlates pos- 
itively with r, which behaviour is opposite to that expected 
in the patchy corona model. We first note that since r <; 2 
in this model, virtually all Compton-reflected photons pass- 
ing through an active region are removed from the observed 
reflected spectrum regardless of the exact value of r. On the 
other hand, a decrease of the observed fi can occur when 
the height of the active region decreases, due to fewer re- 
flected photons being able to escape without hitting the hot 
plasma. This also will soften the spectrum due to the in- 
creased cooling by more blackbody photons emitted by the 



disc intercepted by the hot plasma. The observed opposite 
correlation provides then strong evidence against the patchy 
corona model. 

On the other hand, the data can be well modelled in 
all respects by a geometry with cold clouds inside a hot 
slab (S76; Celotti, Fabian & Rees 1992; Kuncic, Celotti & 
Rees 1997; Collin-Souffrin et al. 1996; Krolik 1998) sur- 
rounded by a cold disc, see Fig. 5d. We assume that the 
clouds are located in the slab midplane and cover a frac- 
tion of f c of the midplane. If there is neither dissipation in 
the cold clouds nor outside seed photons, a solution satisfy- 
ing both the energy balance and spectral constraints corre- 
sponds to f c ~ 0.3, see Fig. 6. The plasma parameters are 
then kT = 51+1 keV and r = 0.95±5J;ig corresponding to the 
half-thickness of the slab (at \ 2 = 45/79 d.o.f.). Compton 
reflection by the cold clouds is attenuated by the hot plasma 
and thus we need the presence of additional Compton reflec- 
tion from an outside matter with Q./2n ~ 0.4 (see Table 2). 
This can occur due to reflection of the radiation of the hot 
flow by an outside cold disc. The covering factor by the cold 
clouds will be f c < 0.3 if additional soft photons from the 
outside cold disc and/or from dissipation in the cold clouds 
(with the luminosity Lbb.intr) enter the hot flow, see Fig. 6. 

As argued by S76, the clouds can be formed by the 
radiation-pressure induced instability of an optically-thick 
disc (Lightman & Eardley 1974; Shakura & Sunyaev 976). 
The clouds are in pressure equilibrium with the hot medium. 

We have also tested a similar geometry with the hot 
plasma forming a central hot sphere (rather than a slab) 
surrounded by a cold disc (e.g. Poutanen, Krolik & Ryde 
1997), see Fig. 5e. If the cold disc does not penetrate into 
the hot cloud, kT = 57± b 3 keV and the radial r = 2.0t°'.l 
(x 2 = 47/79 d.o.f.), and internal dissipation in the cold disc 
is required to provide enough seed soft photons. The model 
also predicts less Compton reflection than observed, which 
problem may be solved by flaring of the outside disc. The 
cold-disc solution of Shakura & Sunyaev (1973) does, in fact, 
correspond to a flared disc, and the illumination by the cen- 
tral hot source will further increase the amount of flaring. 
We note that a central hot sphere is a poor approximation 
to the actual geometry of a h ot accretion flow with sub- 
stantial cooling (see Section 5.3), in which the scaleheight is 
maximized at the outer boundary rather than at the center, 
which effect will also lead to an increase of the amount of 
reflected X-rays. We also find that models with the cold disc 
penetrating into the hot cloud down to Jj 0.7 of the sphere 
radius can also fit the data and satisfy the energy balance, 
but with less or no dissipation in the cold disc. An overlap 
between the cold and hot phases increases the amount of 
Compton reflection in the model and thus reduces the £l/2n 
needed from flaring of the outside cold disc. 

The fl vs. r correlation found by U94 may be explained 
naturally by geometries shown in Figs. 5d, e. In these ge- 
ometries, an increase of the area of the outside soft X-ray 
emitter leads to an increase of the seed-photon flux incident 
on the hot plasma, which in turn leads to more cooling and 
softening of the X-ray spectrum, i.e., an increase of V. The 
same effect leads to more Comton reflection, i.e., an increase 
of Q. 

Summarizing this section, our best models consist of 
a central hot region surrounded by a cold disc (Figs. 5d, 
e). These models require the presence of an additional 
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Figure 6. Relations between the electron temperature, kT, and the half-thickness optical depth, t, of a hot slab, see Fig. 5d and Section 
5.1. The covering factor of cold clouds in the mid-plane of the hot slab is f c — » (upper panel) and f c = 0.3 (lower panel). The solid 
curves correspond to constant ^bb.intr/^hot i i- e -i the ratio of the internally generated and/or external soft luminosity to the dissipation 
rate in the hot disc. These relations do not depend on the amount of pairs. The dashed curves correspond to a constant local hard 
compactness, i (assuming pair-dominated thermal plasma, see Section 5.6). The dotted curves correspond to a constant intrinsic (i.e., 
without Compton reflection) spectral index Y in 4-20 keV range. Crosses are the best fits to the data and elongated confidence contours 
are for 90 per cent confidence for 2 parameters (Ax 2 = 4.61). For f c — > 0, the seed-photon luminosity entering the hot slab is ~ 0.23 
of the heating rate, Lhot- The covering factor of / c = 0.3 is the largest energetically possible consistent with the data. In that case, no 
internal dissipation in the cold clouds and no external soft photons are allowed, and cooling is provided by reprocessed radiation only. 



Compton-refiection component, e.g. from an outer cold disc. 
Among those 2 models, a hot disc with coldclouds insidesur- 
rounded by a cold disc (Fig. 5d) is more likely to be formed 
by an accretion flow with substantial cooling (see Section 



5.3 below). Both models explain the il-T correlation. On 
the other hand, a model consisting of active regions at some 
height above a cold disc (Fig. 5c) is marginally possible, but 
it requires the temperature of blackbody photons lower than 
that inferred from the soft X-ray excess and a truncation of 
the underlying disc, and it predicts an Q-F correlation op- 
posite to that observed. Models with a homogeneous corona 
and with the active regions on the disc surface (Figs. 5a, b) 
are strongly ruled out by the requirement of energy balance. 

5.2 The origin of the soft X-ray excess 

In this Section, we fit thedata in the full 1.2-1000 keV range, 
i.e., including the range < 4 keV sho wing the soft X-ray 
excess. Based on the results of Section 5.1, the most likely 
geometry of the hot plasma is a hot inner disc mixed with 
cold clouds. We first check whether the soft excess can be 
accounted for just by the emission of the cold clouds,which 
we find energetically viable. However, this model fits the 
data much worse than our baseline model of Table 2, A^ 2 = 
+9, even for the maximum possible covering factor of the 
cold clouds (see Fig. 6), / c = 0.3. 

On the other hand, the hot disc is surrounded by a cold 
outer disc beyond a transition radius, R tT (Fig. 5d), which 
contributes to the soft excess. We have thus added to our 
spectral model an integral of blackbody spectra over the 
disc surface at R > R tT . The local color temperature is in 



general larger than the effective one, and the ratio of the 
two for best estimates of the system parameters (Section 
2) is T co i or /T cff ~ 1.8 (Shimura & Takahara 1995). Follow- 
ing Makishima et al. (1986), we neglect here the boundary- 
condition factor J(r) = 1 — (6/r) 1 ^ 2 , where r = Rc 2 /GMx, 
which is a good approximation for a large enough r tr . We 
find that due to the limited energy range of the Ginga data, 
our fits cannot determine the value of Rtr- On the other 
hand, R tI can be obtained from the area covered by the cold 
clouds emitting blackbody spectrum with the fitted color 
temperature, fcTbb- This yields R tT — 6 x 10 7 cm at the best 
fit. The main model parameters are given in Table 3, where 
fcTtr is the cold disc temperature at R tI , and Compton re- 
flection is from the cold disc. The model spectrum is shown 
in Fig. 7. The bolometric luminosity is 4 x 10 37 erg s _1 , 
and the ratio of the blackbody luminosity of the cold disc to 
the total one (but without including the outside reflection 
component) is ~ 0.30. 

We point out that we have neglected the effect of heat- 
ing by the central source on the radial temperature depen- 
dence of the outside cold disc. This may significantly soften 
the power-law part of the disc-blackbody spectrum shown 
in Fig. 7 (e.g. Vrtilek et al. 1990). 

A model alternative to the hot inner disc that was found 



by us marginally possible (Section 5.1) is active regions at 
some height above the disc (Fig. 5c). In this model, the 
soft excess can be accounted for by the blackbody emis- 
sion of the disc region underneath the active region. Due to 
the intensive heating by the X-rays, the local temperature 
can be much above the temperature of a disc with inter- 
nal dissipation only. From the strength of the blackbody 
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Table 3. The main parameters of of our best model, consisting of a hot inner disc with internal cold clouds surrounded by a cold disc, 
fitted to the data set 1. See Figs. 5d, 7 and Section 5.2. Temperatures are in units of keV. 



r kT f c kT hb kTtr X 2 /dof 

0.88±°;Jg 52±J Q.26±°;?t 0.33±g;g 0.21±g;gJ 48.6/81 
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Figure 7. The unabsorbed model spectrum corresponding to the 
fit to the data of our best model, a hot disc surrounded by a 
cold disc, see Fig. 5d. The dot-dashed curve corresponds to the 
disc-blackbody emission, the dashed curve gives the Comptonized 
emission in the hot disc, and the dotted curve gives the sum of 
the blackbody emission of the cold clouds within the hot disc and 
the reflected component from both the clouds and the cold disc. 
The solid curve gives the total spectrum. 



component (Table 2), the characteristic size of the heated 
region of the disc is ~ 10 7 (T co i or /T cff ) 2 cm. However, when 
we approximate the emission of the underlying region of the 
disc as a blackbody, a very poor fit to the data is obtained, 
X 2 = 93/82 d.o.f. This provides one more argument against 
this model. 

Summarizing this section, the most likely origin of the 
soft excess we have found is from an outer cold accretion disc 
with an additional contribution from cold clouds within the 
hot disc, see Figs. 5d and 7. 



5.3 Accretion disc solutions 



5.2 



consists 



The best overall model found in Sections 5.1- 
of a hot plasma slab with the half-thickness corresponding 
to t ^ 1, and kT J> 50 keV. Physically, this geometrical 
model likely corresponds to a hot accretion disc. We note 
that indeed the best-fit parameters of the hot slab are close 
to those predicted by the two-temperature, hot disc model 
of S76. In that model, the local gravitational energy is con- 
verted into the thermal energy of ions (at a temperature 
71 S> T), which is then transferred to electrons by Coulomb 
interactions. The electrons, in turn, radiate away their en- 
ergy by Compton upscattering seed photons irradiating the 
plasma. As found by e.g. Ichimaru (1977), Abramowicz et 
al. (1995) and Narayan & Yi (1995), advection of hot ions to 
the black hole leads to another branch of the hot solution as 
well as it limits the accretion rate possible in the hot flow. 

Detailed properties of the solution in the vicinity of the 
maximum accretion rate are studied by Zdziarski (1998, 
hereafter Z98). That study follows S76 in characterizing 
the flow by a value of the Compton parameter, y = 



49max(r, r 2 ) (e.g. Rybicki & Lightman 1979), which ap- 
proximately determines the X-ray spectral index, F. This is 
equivalent to assuming that flux in the seed photons is such 
that it gives rise to an X-ray spectrum with that T. 

Other parameters of the flow are M and the viscosity 
parameter, a (Shakura & Sunyaev 1973). [The rate of ad- 
vection is characterized by a parameter £ a d v — 1, e.g. Chen, 
Abramowicz & Lasota (1997).] Z98 has obtained that the 
optical depth (corresponding to the scaleheight, H) of the 
flow weakly depends on R, and that r at the local maximum 
rate, M max , and the maximum r (see Fig. 8) are, 



t(M u 



1.22y 3/5 a 2/5 , 



1 c-a 3/5 2/5 

1.54y ' a , 



(6) 



respectively. Those values are reached at r ~ 13, at which 
the rate of dissipation per logarithmic radius interval is also 
close to maximum. At M max , advection carries about half 
the dissipated power and H/R ~ 1. 

The luminosity of GX 339-4 during the observations 
reported is about the maximum observed in the hard state 
(Harmon et al. 1994; Rubin et al. 1998). Thus, we can relate 
the above r with that obtained from fitting the slab model 
in Section ^j]. We assume here an intermediate value of the 
covering factor, f c = 0.1 (see Fig. 6), for which kT = 52^g 
keV, r = 0.82t°;^ ( X 2 = 45/79 d.o.f.), implying y = 0.33. 
Equating the fitted r to r max , we obtain a ~ 1, which is 
relatively large [as it appears to be generally the case in 
luminous black hole systems (Narayan 1996)]. Fig. 8 shows 
the relation between r and rh (= Mc 2 /Le) for those y and 
a; we see that rh max ~ 4. 

The maximum possible luminosity from the hot flow 
corresponds to local accretion rate at the maximum set by 
advection (Fig. 8), for which radial integration implies, 



L hot & 0.15y 3/5 a 7/5 L 



(7) 



(Z98). For the luminosity in Table 2, this relation yields 
Mx <; 2.5M0, which is compatible with our estimates of 
the black hole mass in Section 2. E.g. Lhot — 0.06Le at 
M x = 3M Q . 

Based on our data for GX 339-4 alone, we cannot distin- 
guish which branch of the optically thin solution (Fig. 8) the 
source follows in the hard state at luminosities lower than 
the maximum one. In the a model, the advective branch is 
stable (e.g. Wu 1997) whereas the cooling-dominated branch 
is unstable (Pringle 1976), but the validity of those predic- 
tions has not been tested observationally yet. We point out 
that if the flux variability in the hard state within a factor 
of ~ 50 (see U94) corresponds to the advective branch of 
the disc solution, the average accretion rate in the hot flow 
is within a factor of ~ 2 of M max (compare the left-hand 
dashed and solid curves in Fig. 8). The corresponding large 
(M) requires then the presence of a radiation-driven outflow 
from the companion star for our derived system parameters, 
see Section 2.3. 

The formation of the hot disc might be related to the 
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Figure 8. Accretion disc solutions at r = 13. The optically-thin 
solution (Z98) is given for y = 0.33, and a = 1. We identify the ob- 
served state of GX 339—4 with the maximum of this solution. For 
comparison, we also show the optically-thick solution (Shakura & 
Sunyaev 1973; Abramowicz et al. 1988) computed for the same 
a and M = 3Mq. Solid curves give the accreted fa whereas the 
dashed curves give the part of fa which is locally converted into 
escaping radiation. The remaining part of the locally-dissipated 
energy is advected into the black hole. 



instability of the cold disc in the inner part dominated by 
radiation pressure (Lightman & Eardley 1974; Shakura & 
Sunyaev 1976). If this is the case, the transition radius is 
given by, 



J(r tI ) 



16/21 



60(aM x /M Q ) 2/21 m 16/21 , 



(«) 



(e.g. Svensson & Zdziarski 1994), where a is the viscosity 
parameter in the cold disc, not necessarily the same as that 
in the hot one. For m ~ 4 and a ~ 1, r tr ~ 200. 

Our best-fit model has Ltb/L ~ 0.3, see Section |5.2| 
If the flow is on average ~ 50 per cent advective, L^b/L ~ 
50/r tr (r tr S> 50), which then implies r tr ~ 150, in good 
agreement with equation (^) above. Then, R tT = 6 x 10 7 
cm obtained in Section 5.2 corresponds to Mx ~ 3, in good 
agreement with other estimates of Mx here. 



5.4 The physical state of the reflecting medium 

The region of formation of the Fe Ka line is strongly ionized 
(Section 4, see also U94). At the best-fit ionization parame- 
ters obtained in Table 2, we find that the dominant Fe ions 
are Fexix and Fexx for the first and second observation, 
respectively. Those ions cannot, however, produce the ob- 
served line because of the very strong resonant absorption 
(Matt, Fabian & Ross 1993; Zycki & Czerny 1994). The ob- 
served line is produced by the nonresonant Fe ions ^ xvi. 
We have computed that those ions constitute about 20 per 
cent of all iron in the reflecting medium, which explains the 
weakness of the line. The observed line equivalent width of 
~ 40 eV is then fully consistent with theoretical calculations 
for ~ 3, Q/2tt ~ 0.4 and the obtained ionization state of 
the medium, see Zycki & Czerny (1994), George & Fabian 
(1991) and Matt et al. (1993). That we did not require red- 
shifts of the line center energy nor intrinsic line broadening 



is consistent with the line origin at large distances from the 
black hole, r>l. 

We can also obtain an estimate of the transition radius 
from the ionization state of the reflector. The density, n, 
corresponding to unit Thomson optical depth (where most of 
Compton reflection takes place) can be calculated from the 
vertical structure of the standard cold disc (Milson, Chen & 
Taam 1994). E.g. for Mx = 3M Q and m = 3, n = 1.8 x 10 20 
cm" 3 , 6.0 x 10 19 cm" 3 , 1.8 x 10 19 cm -3 , and 1.9 x 10 17 cm" 3 
at r = 20, 60, 200, and 2000, respectively. If there is a central 
hot disc irradiating the outer cold disc, we can estimate 
the distance between the illuminating source and the disc 
region where most of reflection takes place to roughly equal 
the transition radius. The illuminating luminosity on each 
side of the disc is ~ (l/2)(n/27r)Lhot ■ From that we obtain 
£ ~ 1200, 500, 120, 120 erg cm s" 1 at the 4 values of r above, 
respectively. Thus, we see that r tr i; 10 2 is compatible with 
the fitted values of £ (see Table 2). 

On the other hand, the illuminating luminosity is just 
Lhot in the model of an active region above the disc (Fig. 
5c). From energy balance, we found such an active region 
is l ocated at the height corresponding to its size (Section 
5.1), and this height has to be < R. Then, using the val- 



ues of density from above, we find £ <; 2000 erg cm s~ at 
r Si 100. Thus, the disc ionization expected in the active 
region model is much more than that observed. If there are 
k > 1 active regions with a typical luminosity of Lhot /k, its 
typical size (and thus, the height) must be at least < R/k 1 ^ 2 
to fit on the disc surface, and thus the above estimate of £ 
remains unchanged. This estimate provides one more argu- 
ment against the active region model. We caution, however, 
that the above results have been obtained using our highly 
simplified model of photoionization (see Done et al. 1992). 

Our conclusions here differ from those of U94, who es- 
timated that reflection originates in a region at r ~ 10. The 
difference arises because U94 used the average disc density 
(which is 3> the photosphere density) as well as underesti- 
mated m by using the luminosity in the Ginga energy range 
only and neglecting advection. 



5.5 The role of magnetic fields 

So far, we have neglected in our treatment any magnetic 
fields, which may be present in GX 339-4 (Fabian et al. 
1982; Di Matteo, Celotti & Fabian 1997). If present, ther- 
mal synchrotron radiation will provide additional seed pho- 
tons for thermal Comptonization (Zdziarski 1985, 1986), 
and this process has been suggested to play a major role 
in luminous black-hole sources (e.g. Narayan 1996; Di Mat- 
teo et al. 1997). For kT and r typical of compact objects, 
the synchrotron radiation is self-absorbed up to an energy, 
E s , corresponding to a high harmonic. The emitted spec- 
trum is of the Rayleigh- Jeans form below E 3 , and expo- 
nentially cut off above it (Petrosian 1981; Takahara & Tsu- 
ruta 1982). Zdziarski (1986) has derived an expression for 
E s without taking into account the angular dependence of 
the synchrotron radiation. Here, we modify his expression 
to take into account the effect of angular averaging as cal- 
culated by Mahadevan, Narayan & Yi (1996) to obtain, 



343 a 



In 1 



c 



ln Tn% 



(9) 
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where e B = E s /m c c 2 , e c = B/B cr , the critical magnetic field 
strength is B CI = m 2 c 3 /eh ~ 4.4 x 10 13 G, 



C ~ 70 



rA M exp(l/0) 



3afe c 



2/7 



(10) 



Qf is the fine-structure constant and Am(@, e s /ec) is the low- 
energy and low-temperature correction defined here as the 
ratio of equations (33) to (31) in Mahadevan et al. (1996) 
with coefficients given in their Table 1. The factor C above 
is given in the limit of O -C 1, appropriate for GX 339-4. 

In general, any magnetic field in the source will have the 
strength less than that corresponding to pressure equipar- 
tion. In a hot accretion disc, the la rges t contribution to pres- 
sure comes from hot ions (Section O), and thus B 2 /24tt <> 



nkTi. The ion temperature in the disc is found to be always 
sub-virial, 7] <, 2 x 10 11 K (e.g. S76; Chen et al. 1997). The 
strongest magnetic field is achieved in an inner region, where 
the plasma density is the largest. For estimates, we take a re- 
gion with the size of R ~ 15GMx/c 2 , where dissipation per 
unit logarithmic radius is maximized. In a disc-like geome- 
try, hydrostatic equilibrium implies H/R ~ 1 around Mmax 
(e.g. Z9 8). For r ~ 1 fitted in the slab geometry (Sections 
5.1-5.2), equipartition corresponds then to B ~ 2 x 10 7 G. 
Equations (g)-© then yield E s ~ 10 eV, or e s ~ 2 X 10" 5 . 
{Am ~ 0.1 for the parameters used above.) 

Photons in the resulting self-absorbed spectrum un- 
dergo Compton upscattering in the thermal plasma. The 
spectral index of the Comptonized emission, V ~ 1.75 
(see Table 2) implies certain ratio between the Comp- 
tonized luminosity to that in the seed photons (see Zdziarski 
1986). Assuming the source area of 2irR 2 , the Comptonized- 
synchrotron luminosity is approximately, 



10m o c 3 (2©) 3 " r r+1 
L CS * — ^- r e s 



R 2 



(11) 



where A ~ 2.43 x 10~ 10 cm is the electron Compton wave- 
length. For the source parameters obtained here, Lcs — 
3 x 10 34 erg s _1 . This is 3 orders of magnitude below the 
observed X7 luminosity of ~ 3 x 10 37 erg s _1 . Thus, our con- 
clusion is that the process of thermal synchrotron emission 
is negligible for the formation of the X7 spectrum provided 
the emitting region is similar in size to the region where 
the accretion energy is dissipated. Di Matteo et al. (1997) 
have obtained a larger Lcs in GX 339-4 due to their adop- 
tion of kT ~ 80 keV, a much stronger _B-field than that 
derived above, and neglecting a correction to the rate of the 
synchrotron emission at low kT [as given in Table 1 of Ma- 
hadevan et al. (1996)]. Still, Lcs ~ L for the observations 
considered here would require a source size much larger than 
the region where gravitational energy dissipation is efficient 
and /or magnetic field much above equipartition (equations 

We also note that equations above imply an approxi- 



mate scaling of Lcs/Le oc M- 



(i-r)/2 



Thus, Comptonized 



thermal synchrotron emission from inner hot accretion discs 
around supermassive black holes in Seyfert AGNs (with typ- 
ical T ~ 1.9, 6 ~ 0.2, t ~ 1, e.g., Z97) will be even less 
important than in black-hole binaries. 



5.6 pair production 

If the hot plasma is purely thermal, its parameters of 
kT ~ 50 keV and r ~ 1 imply that e pair production 
is due to photons in the extreme Wien tail of the thermal 
spectrum, and thus it is very inefficient. Indeed, perform- 
ing standard calculations of the pair production rate using 
the model spectrum as in Table 2 for the data set 1 and 
balancing that rate against the pair annihilation rate (e.g. 
Svensson 1984) leads us to the conclusion that there are 
almost no e pairs in the plasma. Specifically, the compact- 
ness, I = LhotfT / Rm e c 3 , required for pair production to be 
able to produce pairs with the fitted r is ~ 10 4 for the hot 
source modelled as a sphere with radius R. This corresponds 
to R ~ 10 s cm, i.e., much less than even the Schwarzschild 
radius of the black hole. 

We have also performed pair-balance calculations for a 
slab, which better approximates the disc geometry. Here, we 
used the model of a sl ab w ith cold clouds in the midplane, 
see Fig. 5d and Section ! 



5.1 



The relevant compactness is then 
the local one corresponding to the power dissipated within 
a cube with size equal to the half-thickness, H, of the slab. 
The results are shown in Fig. 6, where we see that £ ~ 10 3 - 
3 x 10 3 (depending on the covering factor of cold clouds, / c ). 
Approximating the source as a uniformly radiating disc with 
the outer radius R, we obtain R 2 /H = LhotOT /2n£m B (? , 
which for H/R ~ 1 characteristic for hot discs yields 7? ~ 
10 cm, which is the same value as obtained above for the 
spherical geometry, and much less than our estimates of the 
size of the hot plasma. 

Thus, there are no e in the hot plasma in GX 339-4 
provided the plasma is fully thermal. On the other hand, 
it is also possible that nonthermal acceleration of selected 
electrons to relativistic energies operates in the source in 
addition to heating of the plasma. This hybrid model has 
been applied to NGC 4151 by Zdziarski et al. (1996) and 
Johnson et al. (1997) and to Cyg X-l by Poutanen & Coppi 
(1998) and Poutanen (1998). We have applied that model 
to the data set 1. The pair-dominated hybrid model yields 
the same \ 2 as t ne thermal model (fit lc in Table 2) for 
any £ <; 150, and the 90 per cent confidence lower limit 
is £mhi — 70 (for spherical geometry). The fraction of the 
power that is supplied to the accelerated electrons and pairs 
is ~ 0.09 at £ m in, and it decreases with increasing I. This 
^ m in implies R < 10 7 cm, which corresponds to r Si 20 
for Mx = 3Mq. This is in principle compatible with the 
expected size of a hot accretion disc (r ^> 6). 

We have also applied the hybrid model to the data ob- 
tained from the entire OSSE observation of the source, 1991 
September 5-12 (G95), which have much better statistical 
quality than either of our OSSE data sets. We find those data 
are fitted somewhat better, A% 2 = —4, by the hybrid model 
than by the pure thermal one (which gives kT = 62^g 3 
keV at fixed £l/2n = 0.44). The residuals in the thermal 
model show a weak high-energy tail on top of the thermal 
spectrum above ~ 300 keV. This may hint for the presence 
of nonthermal pairs, accounting for the tail. In the hybrid 
model, £ — 150l^J, and the nonthermal fraction at the best 
fit is ~ 0.05. On the other hand, G95 show softening of the 
spectrum of GX 339-4 with time during that observation. 
Therefore, the tail could also be an artefact of fitting the 
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spectrum averaged over a range of plasma parameter by a 
single-component thermal model. 



6 COMPARISON WITH BLACK-HOLE AND 
NEUTRON-STAR SOURCES 

In this section, we compare our X7 spectrum of GX 339- 
4 with those of established black-hole and neutron-star 
sources. Our objective is to determine whether the X7 spec- 
trum of GX 339-4 is indeed similar to those seen in black- 
hole sources but not to those of accreting neutron stars. We 
stress that both black-hole binaries and neutron-star bina- 
ries exhibit two main spectral states, X-ray low (hard) and 
high (soft), and only the corresponding states should be di- 
rectly compared. 



6.1 X7 emission of black- hole sources 

6.1.1 Black-hole binaries in the hard state 

Hard-state X-ray spectra of the archetypical black-hole bi- 
nary Cyg X-l have been fitted by G97. Their fit to the 
average Ginga data of 1991 June 6 with a power-law and 
reflection model (as in Section 4) yields F = 1.591q'q3 an d 
Q/2tv = 0.341qq5 5 which is rather similar to our spectra 
of GX 339-4. X-ray spectra with similar power-law indices 
and moderately weak Compton-reflection components are 
seen from other black-hole binaries in the hard (low) state, 
e.g. Nova Muscae (Ebisawa et al. 1994; Zycki et al. 1998) 
and GS 2023+338 (Zycki et al. 1997). 

The broad-band X7 spectra of Cyg X-l in G97 are also 
similar to that of GX 339-4, and well modelled by a pri- 
mary continuum due to thermal Comptonization in a plasma 
with similar r, but higher electron temperature, kT ~ 100 
keV (and with an additional small spectral component from 
emission of a plasma with t> 1). Soft 7-ray spectra similar 
in shape to that of Cyg X-l have been observed from other 
black-hole binaries in the hard state (e.g. Grove et al. 1998; 
Grebenev, Sunyaev & Pavlinsky 1997). 

Motivated by the similarity, w e te st here the hot accre- 
tion disc model of Z98 (see Section 5.3) against the Cyg X-l 
data. We have refitted the spectrum of G97 correspon ding to 
the peak flux using the hot slab model used in Section 5.1 for 
the primary continuum, and obtained y ~ 0.4 and r ~ 0.5 
(assuming f c = 0.1). Those parameters imply a ~ 0.4 using 
i"(-M max ) of equation (6). The peak luminosity in the hard 
state is about 4 x 10 37 erg s _1 , for which equation (7) im- 
plies Mx iS 10M©. This agrees well with best estimates of 
the mass of Cyg X-l (e.g. van Paradijs & McClintock 1995). 
Thus, it is possible that a hot accretion disc is present in the 
hard state of black-hole binaries, and that its structure de- 
termines both the observed optical depth of the hot plasma 
and the maximum luminosity in the hard state. 
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Figure 9. The ratio of the OSSE spectrum of GX 339-4 from 
1991 September 5-12 to the average OSSE spectrum of the Seyfcrt 
NGC 4151. 



power-law X-ray spectra with F ~ 1.7-2, Compton reflec- 
tion components with Fe Ka lines, and high-energy cutoffs 
above ~ 100 keV (e.g. Z97). The cutoffs can be modelled 
in both cases by Comptonization with r ~ 1 and kT ~ 100 
keV (Z97). 

Here, we point out that the intrinsic spectrum above 
4 keV of GX 339-4 appears virtually identical to that of 
NGC 4151, which, in soft 7-rays, is the brightest and best- 
studied Seyfert AGN. We note here that apart from strong 
absorption, the X7 spectrum of NGC 4151 appears rather 
typical for Seyfert Is, and the OSSE spectra of NGC 4151 
and of all other Seyfert Is observed by OSSE are very 
similar (Zdziarski et al. 1996; Z97). Also, the intrinsic X7 
spectrum of NGC 4151 is very similar to that of a 'typ- 
ical' Seyfert 1, NGC 5548 (Magdziarz et al. 1998). This 
makes models specially designed for NGC 4151 (e.g. Pouta- 
nen et al. 1996) probably unnecessary. The thermal Comp- 
tonization fit (after applying the same correction to solu- 
tions of the Kompaneets equation as in Section 4 above) 
to a joint ROS AT/ Ginga/ OSSE spectrum of NGC 4151 
(Zdziarski et al. 1996) yields V = 1.81±°-°£, kT = 63±\l 
keV , Q/2-k = 0.49t° j6, which are consistent within errors 
with those given in Table 2 for GX 339-4. We find that the 
form of the high-energy cutoff is almost the same in both 
sources, as shown in Fig. 9, which shows the ratio of the 
OSSE spectrum of GX 339-4 (from the entire observation 
of 1991 September 5-12) to the spectrum of NGC 4151 aver- 
aged over all OSSE observations (from Johnson et al. 1997). 



6.1.2 Seyfert Is 

It is of great interest that the X7 spectra of black-hole bina- 
ries in the hard state are similar to the spectra of Seyfert- 1 
AGNs, which similarity in the X-ray regime was pointed out 
by Tanaka (1989). Specifically, both classes of objects show 



6.2 X7 emission of weakly-magnetized accreting 
neutron-stars 

Two types of objects whose X7 emission (if present) is un- 
ambiguously connected with neutron stars are X-ray pulsars 
and isolated neutron stars. Those X7 sources possess strong 
magnetic fields, B <; 10 11 G, and their emission can be easily 
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distinguished from that of black-hole systems, see Finger & 
Prince (1997) and Thompson et al. (1997) for recent reviews. 

On the other hand, weakly-magnetized accreting neu- 
tron stars often show X-ray emission very similar to that 
of black-hole binaries. Some LMXBs classified as black-hole 
candidates have been subsequently found to emit type-1 X- 
ray bursts, leading to their identifications as neutron-star 
binaries. A recent example is discovery of X-ray bursts from 
GS 1826-238 (Bazzano et al. 1997), which source appeared 
as a black- hole candidate in Tanaka (1989). 

Indeed, X-ray emission of X-ray bursters in their hard 
(low) state can be very similar to that of GX 339-4. For ex- 
ample, two Ginga observations of 4U 1608-522 show a power 
law component with V ~ 1.9 accompanied by Compton re- 
flection from a strongly ionized medium with Q/2n ~ 0.5 
in the ~ 2-60 keV range (Yoshida et al. 1993). Similarly, 
GS 1826-238 observed by Ginga has F ~ 1.8 and Compton 
reflection (Strickman et al. 1996). Thus, X-ray spectra alone 
appear not sufficient to distinguish a weakly-magnetized 
neutron star from a black hole. 

On the other hand, typical spectra of X-ray bursters 
observed in the ~ 30-200 keV range by SIGMA telescope on 
board GRANAT and the BATSE detector on board CGRO 
are rather soft, see reviews by Barret & Vedrenne (1994), van 
der Klis (1994), Vargas et al. (1997) and Tavani & Barret 
(1997). In most cases, the spectra can be fitted equally well 
by a power law with (F) ~ 3 or thermal Comptonization 
with (kT) ~ 15-20 keV. There is only one reported OSSE 
detection of an X-ray burster, GS 1826-238, which shows a 
similar V = 3.1 ± 0.5 above 50 keV (Strickman et al. 1996). 
These spectra are all much softer than the corresponding 
spectra of black hole binaries in the hard state [but note 
their similarity to black-hole binary spectra in the soft state 
(Grove et al. 1998)]. 

Unfortunately, there have been almost no observations 
below 30 keV simultaneous with those at J; 30 keV dis- 
cussed above. Therefore, there is no information on the X- 
ray spectral state (high or low) of the bursters during most 
of those observations. However, X-ray bursters are rather 
often found in the low state when observed in X-rays alone, 
and thus we can assume that some of the spectra observed 
at <; 30 keV correspond to the X-ray low, hard, state. Then, 
this implies that the X7 spectra of bursters in the low state 
have spectral breaks at energies significantly below 100 keV, 
or have thermal- Comptonization temperatures of kT Js 30 
keV. This conclusion is indeed confirmed by a ~ 2-250 keV 
observation of 4U 1608-522 in the low state simultaneously 
by Ginga and BATSE, which yields a break at ~ 60 keV 
in the broken power-law model or kT ~ 25 keV, r > 1, in 
the thermal Comptonization model (Zhang et al. 1996). A 
similar spectral break at ~ 60 keV was also seen in a 30-200 
keV observation of 4U 1728-34 (Claret et al. 1994). 

6.3 Spectral black-hole signatures 

As discussed above, X-ray bursters in the low state have X7 
emission with breaks or cutoffs at significantly lower ener- 
gies than both that seen in GX 339-4 and those found in 
black-hole binaries in the hard (low) state in general (Sec- 
tion 6.1.1). When the spectra are fitted by the thermal- 
Comptonization model (with r <; 1), kT ^ 30 keV and 
kT <; 50 keV for neutron stars and black holes, respectively. 



We propose this quantitative criterion to be a black-hole 
signature. 

We see that these two ranges of temperature are not far 
apart, and it is important to test this criterion against avail- 
able data. Those tests should take into account Compton re- 
flection and relativistic effects in thermal Comptonization, 
and use broad-band X7 data. These criteria are not satis- 
fied by thermal-Comptonization fits in Mandrou et al. (1994) 
and Churazov et al. (1994), who find kT = 33 keV and 38 
keV, respectively, in GRS 1758-258, a black-hole candidate 
in the hard state (Sunyaev et al. 1991). Those fits use the 
nonrelativistic model of Sunyaev & Titarchuk (1980), ne- 
glect reflection, and are in a 35-250 keV range only. In fact, 
very similar assumptions lead to a gross underestimate of 
the temperature in Cyg X-l in the hard state, kT — 27 
keV (Sunyaev & Triimper 1979), whereas kT ~ 100 keV 
is obtained from broad-band spectra fitted by a relativistic 
Comptonization model with reflection (e.g. G97). 

As a caveat, we mention that a 40-200 keV SIGMA 
spectrum of a source in Terzan 2 globular cluster appears 
to have no cutoff up to 200 keV and F ~ 1.7 ± 0.5 (Barret 
& Vedrenne 1994; Barret et al. 1991). Such a spectrum may 
imply a Comptonization temperature 100 keV, which is 
characteristic of black-hole sources but not of neutron-star 
ones, whereas the same cluster contains an X-ray burster, 
X1724-308. 

Finally, we point out that the presence of power-law 
emission with F ~ 2.5-3 in the ~ 30-500 keV range, which 
is seen from black hole binaries in the soft state (Grove et 
al. 1998), may not constitute a black-hole signature, in spite 
of a recent claim (Titarchuk 1997). Very similar soft sin- 
gle power-law spectra without detectable high-energy breaks 
are seen from many X-ray bursters (e.g. Barret et al. 1992; 
Goldwurm et al. 1996; Harmon et al. 1996). 



7 CONCLUSIONS 

The results of this work can be divided into two main parts 
as far as dependence on model assumptions is concerned. 
The results of Sections 2, 3, 4 and 6 are relatively model- 
independent. 

In Section 2, we determine the distance to the object 
to be > 3 kpc, with 4 kpc being the most likely value. We 
show the most recent distance determination of ~ 1.3 kpc is 
in error. We also find E B -v — 1.2 and JVh — 6 x 10 21 cm" 2 
as mutually-consistent most likely values. The mass of the 
compact object appears relatively low, Mx ^> 5Mq. 

In Sections 3 and 4, we present our data and show that 
the spectra can be very well fitted by Comptonization in a 
thermal plasma with kT ~ 50 keV and r ~ 1. In addition, 
we find, at very high significance, the presence of Compton 
reflection from an ionized medium and a soft X-ray excess. 

In Section 6, we show that this spectrum is similar 
to that of other black-hole binaries, as well as to those of 
Seyfert AGNs. After comparison with spectra of neutron- 
star sources, we propose that a thermal-Comptonization 
temperature of kT <; 50 keV represents a black-hole sig- 
nature. 

Physical interpretation of our results is given in Section 
5. Here, we study constraints following from the X7 spectra 
and concentrate on accretion disc and corona models. We do 
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not consider, e.g. models with large scale outflows or jets. 
Since the X7 observations were not accompanied by obser- 
vations in any other wavelength, we do not discuss such data 
(e.g. optical, Motch et al. 1983) obtained at other epochs. 
We also do not discuss time variability of the source. 

With these caveats, our best physical model is that of 
a hot accretion disc within ~ 100 gravitational radii sur- 
rounded by a cold outer disc, see Figs. 5d, e. The seed pho- 
tons for thermal Comptonization in the hot disc are supplied 
by cold clouds within the hot disc. The emission of the hot 
disc is Compton-reflected by the outer cold disc. The outer 
disc also emits most of the observed soft X-ray excess. This 
model is in agreement with the spectra, energy balance, and 
ionization balance at the surface of the reflecting outer disc. 
The observed amount of reflection requires that the outer 
disc is flared. 

The hot-disc accretion rate is near the maximum set by 
advection. Based on the spectral fit of the hot slab model, we 
find the viscosity parameter of a ~ 1 and Mx ii 3A4q , which 
mass is in agreement with the dynamical mass determina- 
tion. The hot disc model, which parameters are independent 
of Mx, is also supported by the observed similarity of the 
spectrum of GX 339-4 to those of Seyfert Is. 

We find that e pair production photons in the thermal- 
Comptonization spectrum is negligible and thus the disc is 
most likely made of electrons and ions (although more com- 
plex models with e ± pairs are possible). Also, synchrotron 
emission in the hot disc with equipartition magnetic field is 
negligible as a source of seed photons for Comptonization. 

We can rule out models with a cold disc covered by 
a homogeneous corona or by active regions located on the 
surface of the disc as violating the energy balance. On the 
other hand, the energy balance is satisfied if there are active 
regions at some height above the disc surface. However, this 
model provides poor fits to the spectral data, the predicted 
ionization state of the postulated inner cold disc is much 
higher than that found from the Compton-reflection spectral 
component, as well as it predicts an fi-r correlation opposite 
to that observed. 
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